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Abstract. 

The Causal Dynamical Triangulation model of quantum gravity (CDT) is a proposition to evaluate the path integral over 
space-time geometries using a lattice regularization with a discrete proper time and geometries realized as simplicial mani- 
folds. The model admits a Wick rotation to imaginary time for each space- time configuration. Using computer simulations we 
determined the phase structure of the model and discovered that it predicts a de Sitter phase with a four-dimensional spherical 
semi-classical background geometry. The model has a transfer matrix, relating spatial geometries at adjacent (discrete lattice) 
times. The transfer matrix uniquely determines the theory. We show that the measurements of the scale factor of the (CDT) 
universe are well described by an effective transfer matrix where the matrix elements are labelled only by the scale factor. 
Using computer simulations we determine the effective transfer matrix elements and show how they relate to an effective 
minisuperspace action at all scales. 
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INTRODUCTION 

The aim of the CDT approach is to evaluate the gravita- 
tional quantum amplitude 

G(gi,gf,0:= f ^[^] e'^feMv;^'] (1) 

geometries: gi^gf 

between initial and final geometries {gi,gf}. In this ver- 
sion we do not include matter fields in the theory. We use 
the intuition based on methods of Quantum Field Theory 
to view this amplitude as a path integral over space-time 
geometries. Defining a path integral in this case requires 
solving a number of non- trivial conceptual problems: 

• Definition of time evolution assigns a special role 
to be played by a proper time t for each quantum 
space-time, giving a meaning to the idea of initial 
md final and introducing a time foliation of geome- 
try. 

• At each time t we should define a Hilbert space of 
states - spatial geometries of the Universe. 

• Definition of the measure &[g] is related to the 
choice of the domain of integration (space of ad- 
missible space-times we should include in the path 
integral) and possibly also solving the problem of 
diffeomorphism invariance. 



At the same time one would like to obtain an approach 
which is 

• Background independent - the background geome- 
try may emerge dynamically, but should not be in- 
troduced a priori. 

• Non-perturbative - meaning again that it is not ob- 
tained as a perturbation around some fixed back- 
ground. 

• Has a well-defined infrared limit, described by the 
General Theory of Relativity. 

CDT provides a construction satisfying these require- 
ments [1]. 

• We consider only space-time geometries which ad- 
mit a global time foliation. Causality means that the 
spatial topology of the Universe is fixed during the 
time evolution. 

• We introduce a lattice regularization of geometries, 
assuming that both space and time are discretized. 
The time variable is indexed by an integer time. At 
a fixed time we construct the Hilbert space of states 
defined as a set of states | T) representing triangula- 
tions of a 3d topological sphere using regular sim- 
plices (tetrahedra) with a common edge length a^. 
This definition does not involve coordinates, being 
by construction diffeomorphism invariant. Different 
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FIGURE 1. Building blocks of a simplicial four-manifold 



The parameter A is related to the ratio of the lattice spac- 
ings in time and spatial direction (A = for Us = at). 
For the imaginary time the quantum amplitude has the 
form of a partition function of a statistical ensemble of 
discretized space-time geometries. 

One can observe that the amplitude (partition func- 
tion) can be represented as a matrix product 



G{TiJf-t) = 

^ (7b|^|ri)c(ri)(ri|^|r2). 



(4) 



triangulations |r) correspond to different geome- 
tries, which cannot be mapped onto each other. The 
states |r) satisfy 



{T\T') = ^drr 



(2) 



where Ct is the order of the automorphism group of 
T . In the construction of the Hilbert space we con- 
sider only three-manifolds with the simplest topol- 
ogy of of a three-sphere . This space splits in a 
natural way into a simple sum of spaces labelled by 
the number of tetrahedra N. The number of states 
for a fixed volume N is finite, but large. It grows 
exponentially with N for large N. 
• Tetrahedra at time t are bases of four-simplices 
{1,4} and {4,1} with four vertices at time t and one 
vertex at ^ ± 1 . To form a closed four-dimensional 
manifold we need also four-simplices {3,2} and 
{2,3} with three vertices (triangle) at t and two 
(spanning a link) at r ± 1. Four-simplices have a 
common length of the time link equal to dt see Fig. 
1. The sum (integral) over space-times is regular- 
ized as a sum over simplicial manifolds with topol- 
ogy X [0, 1]. Wick rotation to imaginary time can 
be realized as analytic continuation in dt iat and 
can be performed for each space-time configuration. 

Manifolds can be characterized by a set of global num- 
bers, where and N^^'^"^ denote the numbers of four- 
simplices of a particular type and Nq the number of ver- 
tices. Other numbers of this type can be expressed by this 
triple using topological identities. After a Wick rotation 
each space-time configuration appears in the sum with 

A^o)). where 



the real weight exp(-^(A^|'^'^\A^f '^^ 
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is the Hilbert-Einstein action calculated using the sim- 
plicial structure defined above. Dimensionless coupling 
constants Kq and K4 are related respectively to the in- 
verse gravitational constant and cosmological constant. 



Matrix elements {T\^\T') depend on a number of dis- 
tinct ways to connect geometric spatial states at times t 
andr + 1. 



NUMERICAL RESULTS 

In practice the model cannot be solved analytically, ex- 
cept in the simplest case of 1+1 dimensional space-time 
[2]. For larger dimensionality we are forced to use nu- 
merical methods as a tool to obtain physical information 
about the system. The tools we use are Monte Carlo sim- 
ulations. In the CDT model a very important ingredient 
of the theory becomes the entropy of configurations. For 
a fixed set of Nq^N^^'^^N^^'^"^ the number of space-time 
configurations grows exponentially like 

^(A^o,A^f' ^) - exp«(A^f +A^f '^)). (5) 

The critical parameter depends on Kq and A, which 
means that the parameter A, which originally had a ge- 
ometric origin, plays the role of an independent cou- 
pling constant. The critical parameter renormalizes the 

and 



bare cosmological constant K4 K^^^ = K4 
the model is defined only for K*^^^ > 0. In the limit 
K^^^ 0+ the average total number of simplices goes 
to infinity. This is the limit relevant for the continuum, 
where we may discuss what happens when the discretiza- 
tion effects can be neglected. 

In practice, the numerical approach means that we 
must consider systems with a finite volume. We can re- 
cover the information about the continuum properties 
studying the scaling properties of observables for a se- 
quence of large but finite A^4 ^ 00. This reduces the pa- 
rameter space of the model to a set of two bare coupling 
constants Kq and A. In our simulations we choose peri- 
odic boundary conditions, which on the one hand frees 
us from the necessity to define initial and final geome- 
tries, but on the other hand does not change the physical 
picture, as will become clear below. 




FIGURE 2. Three phases of the four-dimensional CDT. 

Phase structure of CDT 

Using the Monte Carlo program we perform a random 
walk in the space of configurations using 7 elementary 
local moves, which preserve the (local and global) topo- 
logical restrictions on a manifold. The probability to per- 
form a particular move is obtained by the detailed bal- 
ance condition. This is a Markov process with a station- 
ary limiting distribution satisfying 

Prob{^) oc exp(-^(^)) (6) 

Configurations separated by a large number of moves are 
statistically independent and the probability to obtain a 
particular configuration is given by the limiting distribu- 
tion. Expectation values of observables are measured as 
averages in the large but finite set of statistically inde- 
pendent space-time configurations obtained at a particu- 
lar set of parameters Kq and A and (approximately) fixed 

N^'^\ The measurements are repeated for an increasing 

sequence of ' ^ to check the scaling. In the following 
we discuss a very useful observable characterizing each 
space-time configuration. It is the distribution of a three- 
volume N(t) as a function of discrete time t. 

Depending on the position in the {K(),A} plane our 
system appears to be in three physically distinct phases 
[3] characterized by different behaviour of N{t). On the 
Fig. 2 we show a sample distribution of the three- volume 
N{t) for one typical configuration in phases A, B and 
C. Phase A is characterized be a sequence of slices N{t) 
with no correlation between the states at neighbouring 
times. In phase B the time dependence of the distribu- 
tion is squeezed to one time value (one may view it as a 
spontaneous compactification of the time variable). For 
other times the volume is close to minimal. It cannot 
be completely zero, because we choose periodic bound- 
ary conditions and do not allow the volume to vanish at 
any fixed time t. Most interesting from a physical point 
of view is the C phase. The volume profile looks like a 
fluctuation superimposed over a regular classical back- 
ground (the red line on the plot). A typical configuration 
consists of a central hloh and a stalk of cut-off size result- 



ing again from our choice of boundary conditions (peri- 
odicity in time). The red line is the average distribution 
over many configurations with the same volume. We can 
compare distributions for a sequence of volumes (Fig. 3) 
and we find a universal scaling behaviour in the variable 
t = (n — nc) /N^^^^\ with Hausdorff dimension J// ^ 4. 
The plot illustrates the universality of the volume distri- 
bution for the rescaled observable n{t) = N{t) /n\~^^^^ 
plotted vs. t. We expect the averaged scaled distribution 
to be volume-independent. This distribution can be inter- 
preted as a semi-classical limiting distribution of volume. 
Note that in our numerical experiments it is obtained by 
integrating out all other degrees of freedom (details of 
the geometry) except the spatial volume N{t). The anal- 
ysis shows that the averaged geometry scales in a way 
consistent with dimension four. Although this result may 
appear trivial, it is definitely not, since the distribution is 
obtained as the effect of a very delicate balance between 
the entropy of configurations and the physical action. In 
earlier studies, where causality was not imposed, typical 
geometries dominating the quantum amplitude had either 
dn = 2 (branched polymer phase) or J// = oo (collapsed 
phase). 

We can analyze further the properties of the distribu- 
tion and try to fit the limiting curve by an analytic for- 
mula. The effect of the analysis is presented on the Fig. 4, 
together with the fit. The analytic form of the fit suggests 
that the observed geometry can be interpreted as the vol- 
ume dependence inside a four-dimensional ball, in this 
case the variable t plays the role of the azimuthal angle. 
This would indicate that in phase C we see the appear- 
ance of a spherical four-dimensional de Sitter geometry. 

The geometric properties can be analyzed using other 
observables. A useful example is that of the spectral di- 
mension ds. To measure this quantity we analyze the re- 
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FIGURE 3. Scaling of the volume profiles. 



Measured data 




FIGURE 4. The averaged volume profile and the analytic 
formula used in the fit. 



turn probability in the diffusion process on the geometry, 
as a function of the diffusion time a [4] . If the geometry 
was regular we would expect 



1 



ds 



/logP(cj) 
Jlogcr 



(7) 



with a constant ds = dn- The figure shows the observed 
behaviour of ds obtained by averaging over many starting 
points of the diffusion process and over many configura- 
tions. The plot on Fig. 5 shows that ds is not a constant, 
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FIGURE 5. Scale dependence of ds. 

but depends on a suggesting a scale dependence of the 
effective geometry, ranging between two at short scales 
and four at large scales. This illustrates the quantum char- 
acter of geometry. A similar property was discovered in 
other approaches to quantum gravity (c.f. e.g. [5]). 

The geometry presented above was measured at a 
particular point on the {ko,A} plane. When the values 



of the bare coupling constants are changed inside phase 
C, the qualitative behaviour remains the same, up to a 
finite change in the scale. We are particularly interested 
in the critical behaviour near phase transitions. 

The qualitative behaviour of the phase structure of 
CDT was found to have strong similarity to the phase 
structure predicted by Hofava-Lifschitz gravity ([6] and 
[7]). The simplest check of the analogy was to measure 
the order of the phase transitions, which we found to be 
first-order for the A-C transition and second order for the 
B-C transition [8]. 



Effective action in de Sitter phase 

The regular semi-classical distribution of spatial vol- 
ume observed in phase C (de Sitter phase) suggests that 
it reproduces a saddle point of some effective action of 
the spatial volume V{t) or, equivalently, of the scale fac- 
tor a{t) = V{ty^^. A natural candidate for such an action 
is the mini- super space action 



S = 



1 



247rG 



J dt ^/git(^ 



V{t) 



(8) 

where t takes the continuum value and A plays the role 
of Lagrange multiplier, necessary to fix the total volume 
to some target volume V 



v = Jdt v^y(0. 



(9) 



Here sets the scale in the time direction. In a discrete 
setup we may expect this action to take the discretized 
form 



1 



(10) 



1/3 



although it could have a more complicated form. Here 
Nt=N{t). 

Measuring the covariance matrix of volume fluctua- 
tions around the semi-classical distribution and invert- 
ing this matrix we can determine the matrix of second 
derivatives of the effective action, assuming that higher- 
order terms (higher than second order in fluctuations) can 
be neglected. This method was successfully applied in 
[9]. Indeed the form (10) of the effective action was con- 
firmed, at least in the range of large volumes, permitting 
us to determine the physical parameters F and fl. The pa- 
rameter /i was particularly difficult to measure, since the 



corresponding term in the action (after differentiating it 
twice) falls off very fast with the spatial volume. On the 
other hand, small volumes can be expected to be (and 
are in fact) very sensitive to finite- size effects and lattice 
artefacts. 

The numerical experiment described above produces 
values of the physical parameters {r,/i} as functions of 
the bare couplings {/Co,A}. These physical parameters 
have a direct interpretation in terms of the gravitational 
constant G (up to the dimensionful parameters at and Gs). 
From a practical point of view the determination of these 
parameters becomes more difficult near the phase tran- 
sitions, where we observe a critical slowing-down and 
large finite- size effects. As we observed, the neighbour- 
hood of the phase transitions is particularly interesting 
from a physical point of view. Approaching these lines 
we see that the parameter l/F ^ (or F ^ oo) [10], 
which can be interpreted as the limit where the lattice 
spacing approaches zero. In this limit we may hope to 
observe genuine quantum effects of gravity. 

TRANSFER MATRIX 

The form of the effective action (10) suggests a formal 
decomposition of the quantum amplitude (4) in the sim- 
plified form 

G{No,Nt)= (11) 
^ {No\^\Ni){Ni\^\N2) {Nt-i\^\Nt) 

where we have introduced the effective projection opera- 
tors 

p{N) ^ \N){N\= ^ \T)C{T){Tl (12) 

p{N)p{N^) = Wp(A^) 

on the space of states with a fixed volume N. The projec- 
tion operators behave as genuine projection operators on 
a single state \N). They can be used to study the proper- 
ties of the CDT geom etry, a ssuming that t he eleme nts of 
the transfer matrix ^/C^ {Ti\ ^ \Ti^ i) ^/^fi^ in (4) 
depend only on volume. We can check to what extent this 
is true. 

In the proposed approach [11] we determine directly 
the elements {Ni\^\Ni^i) using numerical simulations 
of periodic systems with very small time extent. The 
method is based on the observation that terms in the sum 
(11) have the interpretation (up to a normalization) of the 
probability to measure a particular sequence of volumes. 
For a system with periodicity 2 and periodic boundary 
conditions 



By measuring the number of times a particular set 
{A^i , appears in the simulation we determine the ma- 
trix element {Ni \^\N2).ln practice the method is more 
complicated, because we also want to study a particular 
range of {A^i , A^2}- Details of the method are explained in 
[11]. On the Fig. 6 we show the logarithm of the transfer 




FIGURE 6. The logarithm of the transfer matrix elements 
plotted as a function of Ni and N2 for a sequence of ranges. 

matrix, obtained by gluing together the results of mea- 
surements at the neighbouring ranges of volume. On the 
plot we see the Gaussian behaviour of the off-diagonal 
kinetic term and the diagonal potential terms. The kinetic 
term corresponds to the first line in (10) and the poten- 
tial term to the second line. We can easily measure the 
parameters of the effective action and find consistency 
with the values determined by the indirect method de- 
scribed before. The advantage of the new approach is 
a much smaller numerical error and at the same time a 
much shorter computer time needed to perform the mea- 
surements. 

The presented plot corresponds to one particular point 
on the {/Co, A} plane, well inside the de Sitter phase. We 
are currently measuring the behaviour of physical param- 
eters {F, 11} in the whole range of the C (de Sitter) phase, 
in particular near the phase transitions. Preliminary re- 
sults confirm that 1 /F ^ at the A-C transition and in- 
side the A phase and indicate that 1 /F changes sign at the 
B-C transition lines. Details of this behaviour are crucial 
to determine and understand the critical behaviour near 
the phase transition and the critical scaling properties of 
the model. Particularly interesting is the perspective to 
study the neighbourhood of the triple point, where the 
three phases meet. 



^{NuN2) oc {Ni\^\N2){N2\^\Ni) = {{Ni\^\N2))^ 

(13) 



CONCLUSIONS 

The CDT model allows us to study properties of the lat- 
tice regularized quantum theory of geometry in a Wick- 
rotated formulation (imaginary time). Obvious questions 
about the full properties of the model under analytic con- 
tinuation to real time remain open. Some features of the 
model are however common to a formalism with real and 
imaginary time. One important property is the crucial 
role played by the entropy of configurations, an aspect 
which is usually not appreciated in mini- super space-type 
models. Our approach is based on integrating out all de- 
grees of freedom, apart from a finite set and permits us 
to study the true effective model of the scale factor. 

The new method to analyze the properties of the 
model, discussed in this article, gives rise to the hope 
that approaching the boundaries of de Sitter phase will 
be possible. 
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